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Development of a general time-dependent absorbing potential for the CATM 
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The Constrained Adiabatic Trajectory Method (CATM) allows us to compute solutions of the time-dependent 
Schrodinger equation using the Floquet formalism and Fourier decomposition, using matrix manipulation 
within a non-orthogonal basis set, provided that suitable constraints can be applied to the initial conditions 
for the Floquet eigenstate. A general form is derived for the inherent absorbing potential, which can reproduce 
any dispersed boundary conditions. This new artificial potential acting over an additional time interval 
transforms any wavefunction into a desired state, with an error involving exponentially decreasing factors. 
Thus a CATM propagation can be separated into several steps to limit the size of the required Fourier basis. 
This approach is illustrated by some calculations for the molecular ion illuminated by a laser pulse. 



I. INTRODUCTION 

The development of new tools for solving quan- 
tum dynamical problems remains a very active field. 
One relatively modern method, the Constrained Adia- 
batic Trajectory Method (CATM), involves a Floquet 
processing^— of the time-dependent Schrodinger equa- 
tion (TDSE) , using some form of artificial absorbing po- 
tential to connect a single Floquet eigenstate to the wave- 
function of the system (namely the constrained Floquet 
state). The Floquet formalism is frequently used to de- 
scribe interactions between molecules and fields because 
it makes it possible to separate the quasi-periodic oscilla- 
tions of the electromagnetic field from the slow, adiabatic 
variations of the envelope, allowing an efficient descrip- 
tion of singular phenomena such as exceptional points 4 . 
In this paper we choose a slightly different approach, in- 
volving a rigorous periodicity, because all the time inter- 
val to be treated is included in the fundamental period, as 
in the (t,t') approach*. The dynamical problem is thus 
transformed into a "static" problem, in the sense that 
the time will be included in an extended Hilbert space 
which is the product of the usual Hilbert space with the 
space of T-periodic functions. 

This approach is different from the traditional one, 
in which the Hamiltonian is considered as constant over 
small time steps, with propagation of the wavefunction 
using a sequence of differential approximations. In the 
CATM the dynamical integration is completely replaced 
by the search for one eigenvector of a large matrix. 
That can be a difficult task, but many methods exist to 
solve this problem. We may use in particular the time- 
dependent wave operator theory^i, since the absorbing 
potential dilates the Floquet spectrum, isolating the re- 
quired Floquet eigenvalue in a favourable manner—. 

We briefly summarize the main ideas of the CATM, 
already detailed in two articles^. If the Hamiltonian of 
the unperturbed system is Ho(q), with eigenvalues and 
eigenstates {Ej , \ j)}, and W(q,t) is the time dependent 
part of the Hamiltonian, then we can work with the Flo- 



quet Hamiltonian in the extended Hilbert space: 

H F (q, t) = H (q) + W(q, t) - ih^- + V(t). (1) 

at 

Indeed, eigenvectors of Hp are linked directly to solutions 
of the TDSE. The fundamental period T is here chosen as 
the entire duration [0, To] of the physical interaction W(t) 
plus an artificial time interval [To, T] (T = Tq + AT), on 
which the absorbing potential V(t) will act. 

The Floquet states are indexed with a double index^ 
linked to the molecular eigenstates (from Hq, j «-» |j)) 
and to an FBR Fourier basis set (n <H- (t\n) = e 
These eigenstates {|Aj, n (<7,i))} are defined by 

H F \X j , n {q,t))=E Xj jX j , n (q,t)). 



-mint T 



(2) 



If one of these eigenstates projected at t — is pro- 
portional to the selected initial value of the wavefunc- 
tion \l/(0), then the knowledge of this eigenvector im- 
plies automatically that of the wavefunction driven by 
the Schrodinger equation. This follows from the rigorous 
expansion^: 

|*(g,t)> = X>i,n=o(<?,t = 0)1*0?, t = 0)) 

3 

x e -iE, j>n=0 t/H lXjfiM (3) 

which can in this case be reduced to one single term: 



l*M)> 



X(q,t)) 



with 



\V(t = 0)} = a(t = Q|A) 



(4) 



(5) 



a ' Electronic mail: Arnaud.Leclerc@utinam.cnrs.fr 



where the overlap a — (X(t = 0)\^f(t = 0)) is a complex 
number. The time-dependent absorbing potential added 
to the real Hamiltonian ensures the approximate validity 
ofEq.©. 

In previous articles, two different forms for the time- 
dependent absorbing potential were proposed. With this 
artificial potential the CATM was able to deal with many 
simple examples. For a 2-level system, a non-diagonal 
matrix able to reproduce any initial condition has al- 
ready been presented^. Except for this particular case, 
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until now the method could only integrate the dynamical 
evolution of a wavefunction initially projected onto one 
channel j = I of the non-perturbed system 1^(0)) = 
because the absorbing potential was constructed to ab- 
sorb all channels except the initial one. It was repre- 
sented by a diagonal complex matrix, acting over the 
additional time interval [Tq,T] and taking the form 



v(t)=j2 v o P t(t)\j)(j\ 



(6) 



with the three implicit conditions that V^ pt (t) is strictly 
null over [0,T ], Ss(V£ pt (t)) is negative over [T ,T] and 

| $(V c pt (t))dt\ > |9f(£ A )|(T - T ). This form of 
potential allows us to propagate an initial condition 
1^(0)) = \l) since, owing to the periodic properties of Flo- 
quet eigenvectors, the final condition |A(T + AT)) oc \l), 
obtained with the introduction of V(t), induces the initial 
one |A(0)> oc |Z). 

The subject of this paper is to find a general formula 
for an absorbing potential able to constrain the Floquet 
eigenvector to any required boundary condition. With 
this improvement the CATM is able to treat dispersed 
initial wavefunctions and can also propagate the solu- 
tion over very long time intervals (long laser pulse, train 
of pulses...), because the time interval can now be cut 

into smaller pieces [0, Ti], [Ti, T 2 ], So we need a 

new absorbing operator, acting in the extended Hilbert 
space and capable to transform any quantum state into a 
given different state, under the evolution dictated by the 
TDSE (for instance, able to transform \^(t = T 2 )) into 
|*(T 2 + AT)> = |*(Ti)». 

To this end, the paper is organized as follow. In section 
HTl we construct a treatment based on a non-orthogonal 
basis set, as a "memory" of the initial state to be re- 
produced. Unfortunately this intermediate result is not 
exact and an improvement with a correction term is pro- 
posed in section Hill After that the asymptotic form of 
the Floquet eigenvector becomes conclusive. Some tests 
are made in section IIVI on submitted to an intense 
laser pulse. For this same system we compare two prop- 
agation schemes : a one-step scheme with the old version 
of the absorbing potential and a multi-step scheme with 
the new form of the absorbing potential. 



II. A HEURISTIC DEFINITION OF THE ABSORBING 
POTENTIAL 



It is difficult to find directly an absorbing potential 
which acts selectively on each channel to impose the de- 
sired boundary conditions A(T + AT) = ^(0). However, 
we can make profit from the previous definition of a V 
absorbing on all channels except one, by taking the ini- 
tial wavefunction \ty(t = 0)) in place of the state \l), at 
the expense of constructing a non-orthogonal basis set. 



A. Non-orthogonal basis set 

A non-orthogonal basis set for the Hilbert space % is 
defined as 



li) = 



li> W 

|#(0)) for j = l 



(7) 



where j is a global index for the free molecular eigen- 
basis. I is chosen among the {j} so that |(i|\lj(0))| is a 
maximum, in order to avoid numerical difficulties (such 
as divisions by too small numbers, because the later re- 
sults will include (/^(O)) -1 ). Thus the non-orthogonal 
basis is identical to the primitive orthogonal basis, except 
for one single vector \l) which is replaced by the initial 
desired wavefunction 'I'(O)). In a sense, the initial wave- 
function is "kept in memory" as one of the basis vectors 
by using this definition. The new basis for the extended 
Hilbert space K = T-L® ^(T) is thus: 



\U) 



(8) 



where \ti) represents the DVR grid basis on time (the 
small ti being collocation points). This can be written 
equivalently by Fourier transform on the FBR basis \ j) £g> 
|n). 

The transformation matrix B from the orthogonal ba- 
sis set {|fc) <g> \ti)} to the non-orthogonal one {\j) ® |i m )}, 
such that \j) ® I**) = ^2 km B k jS mi \k) ® \t m ), is consti- 
tuted by a series of identical matrices, distributed along 
the main time diagonal S m i- Inside one block ti it is 



B kj = 



S k j if j + I 
<fc|*(0)) if j = l. 



(9) 



For the inverse transformation matrix, one obtains: 



Skj if 3 I 



W*(Q)> 
WW) 



if j = I and k ^ / 



(10) 



-mm 



if j = I and k = I. 



These transformation matrices will be useful in the next 
sections. 



B. First conjecture 

Our first approach is to use the expression (j6]) with- 
out modification within the non-orthogonal basis set, so 
that an absorbing potential is placed on each channel 
\j) ^ \l) but the channel \l) = |*(0)) is kept intact. We 
presume that this definition will produce the same ef- 
fects as previously, by giving an asymptotic absorption 
of components on all channels except the one which is 
proportional to the desired initial wavefunction. Such a 
potential expressed in the non-orthogonal basis of Eq.([7]) 
is (see Tab|T|: 



v = EE li>l^>M3l K c pt (ti) 



(11) 
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TABLE I. Structure of the matrix which represents one block 
ti of the absorbing potential defined in Eq.(QTJ within the non- 
orthogonal basis. This matrix multiplies the factor V p t (ti) in 
Eq. ljlip . All channels are expected to be absorbed except the 
one (I) "containing" |*(0)>. 
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(column I) 



Ol*(0)> 

~ VWW) 
(row I) 

(j|*(0)> 

<i|*(0)> 
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TABLE II. Structure of the matrix which represents one block 
ti of the absorbing potential in the orthogonal basis. This 
matrix multiplies the factor V£ pt (ti) in Eq. (|13|) . 



where V^ pt (t) is different from zero only during the added 
artificial time interval t £ [To,T], 

Further use of this non-orthogonal basis would create 
some complications in the Floquet Hamiltonian repre- 
sentation owing to the non-diagonal expression for Hq, 
so we prefer to come back to the orthogonal basis for 
practical application. The potential must thus undergo 
the following transformation: 



This is no longer a diagonal matrix, but now is, for each 
I 



\non—orth.b. r> — l 



U block (see TabHTJ): 
(ti\V\t i )=V c pt (t i ) x 



EE 



(1 - Sji) S jk + S k i x 



<J|*(0)> 



(13) 



If this definition is appropriate the Floquet eigenvector 
calculated with the modified Hamiltonian should satisfy 
the boundary conditions. In the most general form: 

(j\X(t = 0)) = a x <j|*(0)) Vj (14) 

where a is a complex proportionality factor. 

C. Consequences of the first definition 

It is now necessary to analyse the effect of such a poten- 
tial on Floquet eigenvectors, during the additional time. 
We solve analytically the first order differential equation 
which drives the Floquet eigenvector components (j\X(t)) 
over the time interval [To , T] where the absorbing poten- 
tial is different from zero, while simultaneously all cou- 
plings are reduced to zero. Remembering that Floquet 
eigenvectors satisfy, in the extended Hilbert space /C, the 
equation 

(H F - E X )\X) = 0, (15) 
projecting on the time and writing (t|A) = |A(t)), we get 

^-\X(t)) = ±[H + V(t)-E x }\X(t)). (16) 

For the time being, V is the operator defined in Eq. (ITT1) . 
Ea. (lT6|) with V(t) expressed in the orthogonal basis (see 
Eq. (fT5jl ) has an analytical integral solution involving E\ 
as an undetermined parameter (implicit solution). 

Using the notations = Xj and | \T/ (0) ) = tp®, one 
can write : 



d_ 

Iff 



1 



\ 



( 1 (col.l) 



Ex + V % t {t) x 



1 



(lin.l) 



1 



Xj(t) 



(17) 



Thus for all t within [To,T] we have 
1 



d_ 

Iff 



^(E, - E x )X,(t) 



1 



which directly gives 

In the particular case of the line j = I, this equation 



I 

becomes simply : 
8 



(19) 



A(t), = X t (t = T )exp ( -{E x - Ei)(t ~ T ) j , (20) 
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while Eq. (fT8)) for j ^ I is somewhat more complicated. 
After some elementary simplifications one obtains a for- 
mula describing the Floquet eigenvectors behaviour when 
t€ [T Q ,T] for j^l: 

Vj ^ / A 3 -(t) = Xj(T ) x A,^-^- y ^(*'» dt ' 



To 



x e * St b («* - ^ - v?, t (f ) ) dt' _ ^ 

(Note that Eg . ([20)1 can also be obtained as a particu- 
lar case of Eg . (|2"Tj) ~) . At this point we can ask whether 
Eo. (l2~T1) is consistent with the equality |A(i — T)) = 
a |W(0)). One can first transform the expected propor- 
tionality relationships ([5]) into an equality: 



|A(t = 0)) = |A(i = T)) 
by introducing the vector 

\X(t)) = \X(t)) x 



l*(0)> 



Mi) 

with the components (see Eg . (f2Tj) ) 

'A(t)+^° x B{t) ifj^l 



Xj(t) = 



0? if 3 = 1 



(22) 



(23) 



(24) 



and with 



M{J-o) 



(E,-Ej-V°(t'))dt' 



(25) 



B(t) 



i 

xe 



T 



(26) 



Thus, the eigenvector will satisfy the expected boundary 
conditions (|22|) only if we have 



(27) 



J A(i = T) ~ 
1 £(i = T) ~ 1. 

The first condition implies that J^T 3f(-E; — Ej — 
Vg pt {t'))dt' is large, creating an strong exponential de- 
crease. However a detailed analysis of the second func- 
tion B remains difficult, because of the double interlocked 
integral. Here it is impossible to indicate clearly when the 
second condition is respected. However, numerical tests 
reveal that there is a non-perfect control of the boundary 
conditions. It thus seems necessary to introduce a cor- 
rection term to the definition of V(t) in order to enforce 
the required boundary conditions. 



III. A RIGOROUS DEFINITION ABLE TO REPRODUCE 
ANY INITIAL CONDITION 

To solve this problem rigorously, we follow an inverse 
solution. Considering the results of section IIIC[ it is 



possible to change the solution \X(t)) slightly to obtain a 
proportionality to the initial conditions, and to find after- 
wards what should be added to the absorbing potential so 
as to obtain this modified result. That is, we start from 
the solution to go back to the appropriate equation. 



A. Ideal solution 

If we add in Eq. ([2T]l a factor (V^ pt (t) + Ej - .Ej) at the 
"feet" of the integral, instead of V° pt (t), i.e. 



(28) 



& a {Ej - S / +y o c pt (i'))e i/ ^ (B ^ £l+V/ - t(t " ))£i *"di'(29) 
then this integral simplifies in the following way: 

lT ( V o P t(t') + Ej - E l )e* S i {Ei - El+v °* t{t '' ))dt " dt' 
= f 4 i_d_ A lT (Ej-Ei+v^ t (t"))dt" 

J T i dt' 

and the solution for Xj (t) takes the form: 

Xj{t) = Xj(T )e*ti ( E *- E j-v o c pt {t'))dt' (31) 

As in the previous section, we introduce 



\X(t)} = |A(t)> x 



(32) 



with components 



where 



A; (To) 



D(A = 1 _ e i P To (E l ^E ] ~V^ t (t , ))dt^ 



(35) 



The asymptotic condition |A(T)) = 1^(0)) becomes eas- 
ier to satisfy: we require 



C(t = T) ~ 
D(t = T) ~ 1 



(36) 



this is equivalent to 

eif^'-^-WW ^0 (37a) 

«• / 9f(£,-B J --V4 t (t))dt»l. (37b) 

J To 

At this stage, these corrections (Eqs. |2"8]|2T))) are still 
inconsistent with the selected absorbing potential. We 
see now what change we should make in the absorbing 
potential so as to produce a favourable conclusion after 
integration of the differential equation (TTI 
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B. A correction term to the first definition 

Progressively tracing back the solution, we deduce that 
all we have to do is to introduce a corrective energy term 
(Ej — Ei) inside the absorbing potential. The matrix 
representation is shown on Tab. IIIII for each Floquet 
block £j. The definition (fl"3"| is thus replaced by 



(*<|v|ti> = £E 



(l-fy) [ $jkV c pt (U) + 6 kl x 



-UMQ)) 



(V c pt (t t ) + E 3 - E t ) 



\j)(k\- 



(38) 



V % t (ti) 










(column / 










v: pt (n) 





<3l*(0)> 
<1|*(0)> 




+ Ej - Ei) 




























{) 




(row I) 
















<iI*(o)> 

(1|*(0)> 




+ Ej - Ei) 





TABLE III. Matrix representation of one block U of the ab- 
sorbing potential (second version) within the orthogonal basis 
set. 



Eq. (|16)) for the Floquet vector components is replaced 



by 



%(t) = ^(Ej-Ex + V^ityXj® 



dt 



i h 



L(V c pt (t) + Ej - E t ) 



kit). (39) 



After integration, the formulae can be simplified, as in- 
dicated in Eq. (f30|) . Finally we obtain the ideal solution 
of Eq. ([T4l) . solving the differential equations (f39j) in pres- 
ence of the correcting term and taking t = T: 

Xj(T) = A 3 (T )e- |/ ^ (£A - B ^ y ^ (t ' ))dt ' 



-A 



(To)f (e* 



(E x -Ei)(T-T ) 



(40) 



The desired asymptotic behaviour 
(j\X(t = 0)) = ax 



with a = 



A; (To) 



HEi-EOiT-To) 



(41) 



is obtained if the factor e K ^ {El E ' v °?* {t ' ))dt ' tends 
properly to zero. It means that the integrated surface of 
V£ pt (t) over the time dimension must possess a large neg- 
ative imaginary part, so that a strongly decreasing real 
exponential appears, becoming almost zero when t = T: 

(42) 

In this formulation the CATM will also run with any 
initial superposition of states. 



IV. TESTS ON H£ SUBMITTED TO AN INTENSE 
LASER PULSE 

A. Illustration of the integration scheme 

Using this new form for the absorbing potential, we 
have deal with a test-system. Many propagations were 
made on the example of the molecular ion, mod- 
elled by its first two electronic surfaces 2 E+ and 2 S+ 
and illuminated by an intense laser pulse of total du- 
ration 750a. u. (18fs) with a carrier wave frequency 
0.335a.u. (i.e. a wavelenght of 136nm) and an inten- 
sity of 10 13 W.cm~ 2 . This pulse is represented in FigfT] 
The initial vibrational state was chosen as v = 2 (second 
excited state) because beginning in this state with this 
precise carrier frequency and this relatively high inten- 
sity leads to an electromagnetic trapping—. This choice 
makes the test more significant since the CATM algo- 
rithm should recur strongly non-linear phenomena. 




FIG. 1. Pulse with a total duration of 750a.u., carrier wave 
frequency 0.335a.u. and intensity W 13 W.cm~ 2 . 

Indeed, FigJ5] shows that in these intense pulse condi- 
tion the dissociation remains small (final Pdiss — 7.14 x 
10~ 2 ). The initial state keeps the main part of the pop- 
ulation, while states 1 reaches 5% at its maximum. The 
fundamental state population shows important oscilla- 
tions between 5 x 10 and 10~ 2 . Other bound states, 
with higher energies than the initial state, show many os- 
cillations in population, with a frequency which increases 
gradually as one moves up in energy. 

Since our principal aim is to test the CATM method 
in its more general formulation, we divide the time inter- 
val into N s steps and we realize a CATM propagation for 
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FIG. 2. Dissociation (black) and transition (colours) proba- 
bilities P{\v)) = |H*(t))| 2 for H% submitted to the pulse of 
Figffl The initial state was the second exited state v = 2. 



each step [0, Ti], [Ti,T 2 ], [T 2 , T 3 ] . . . using the general for- 
mula of Eq. (j55|) . Each interval [Ti,T i+ i] is successively 
considered as the physical interval to be submitted to 
the CATM procedure, which requires for each step the 
addition of an artificial time interval of duration AT, 
on which the absorbing potential is introduced. We just 
have to make the correspondance t' = t — Ti (i.e. T ■<-» 0, 
T i+ i o T m a U formulae in section Hll| . 

The first step is a propagation issuing from a sin- 
gle eigenstate |^(0)) = \v = 2), during which the wave- 
function is dispersed over all the basis set. Then for 
the following steps, the CATM must ensure the cor- 
rect reconstruction of the totally dispersed initial states 
(|$(Tj)) = J2 v c v\ v ))- These more general problems are 
rigorously solved provided that Eq (|37bj) is satisfied. This 
is the case if we use the modified expression for the 
absorbing potential; however the supplementary terms 
(Ej —Ei) that it includes create numerical difficulties. 

For each step, over the additional time interval AT, 



terms 



{;|^(o)) (Ej — Ei) are present in the column num- 
ber I of the Hamiltonian, where ^(0) is the wave func- 
tion at the initial instant of the step. Given that the 
propagation is made using numerous Fast Fourier Trans- 
forms (FFT) , these terms must be multiplied by an Heav- 
iside function of time which is zero everywhere except 
on the additional interval of duration AT. We thus en- 
counter some numerical difficulties because the Heaviside 
function is discontinuous and the FFT needs very high 
basis frequencies to describe correctly such discontinu- 
ities (Gibbs phenomenon). Moreover, the electric field 
must be zero over the artificial time interval, so as to 
avoid transitions during the absorption phase. To limit 
this problem we choose the time step separation care- 
fully, corresponding to nodes of the electric field. This 
is the easy way but it does not restrict the generality of 
the approach. In principle, it would be also possible to 
choose other step positions with non-zero values of the 
electric field, but in this case we would have to ensure 



the continuity with an artificial return to zero during the 
additional time interval. 

In this framework, the determination of one Floquet 
eigenvector was realized with a wave operator algorithm 7 , 
starting from a test vector which is constant over the time 
step and modifying it iteratively with a RDWA (Recur- 
sive Distorted Wave Approximation) procedure. 

To illustrate the principles of the method, Fig. [3]shows 
results for a CATM calculation made with an 8-step prop- 
agation. At each step the absorbing potential constrains 
the wavefunction components to match the values ob- 
tained at the end of the previous step. The small sup- 
plementary pieces of curve correspond to additional time 
intervals during which the absorbing potential is present. 
FigHJis a "zoom in" on these results for one given step. 
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FIG. 3. Evolution of some components of the wavefunction on 
bound states of , |(t>[\P(t))|, v = 0, 1, 3. Illustration of the 
CATM scheme with a decomposition of the propagation into 
8 big steps, ~ 94a. u. are treated for each CATM step. Each 
point in the figure corresponds to a grid point of the time- 
discretization. Grid points obtained during additional time 
intervals are kept here to help in understanding the diagram. 



B. The stability of multi-step propagations 

We now compare three different calculations made 
using the following parameters: run A, one step with 
N = 1024 Fourier basis functions; run B, two steps with 
N = 512 and run C, four steps with N = 256. For 
each calculation the absorbing potential amplitude was 
chosen sufficiently large to correctly reproduce the initial 
conditions, so that the jumps between each step and the 
differences between the A, B, C curves become both sta- 
ble and negligible. In other words, the time integral of 
the absorbing potential was always sufficient to stabilize 
the numerical values of the calculated probabilities in the 
three different configurations. 

Previous calculations made on this system using a one 
step procedure^ have proved that the CATM results are 
very precise when compared with ones obtained by using 
standard wave-packet propagation techniques. Thus run 
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FIG. 4. "Zoom in" on the fourth step. Evolution of some 
components of the wavefunction |(v|W(t))[, v = . . . 14 (same 
example as in the previous figure). The real interaction is 
extended with an artificial interval over which the required 
initial conditions are progressively recovered. 



FIG. 6. Evolution of the relative difference on the transition 
probability \(v = 13|W(t)}|, between runs A (1 step), B (2 
steps) and C (3 steps), | g ~ | (green triangles) and | C ~ A | 
(red points). 

0.01 r— ~T ■ — i ■ — i ■ — T ■ — I ■ — , ■ — 1 



A is chosen as the reference and relative differences are 
successively computed between A and B and between 
A and C. Their evolutions are shown in FigJS] Since 
the time when the dissociation probability is no longer 
negligible (t ~ 125a. u.), differences become rapidly small 
and stay so. We cannot see any significant increasing at 
the moment of step changes (t = 187. 5a. u., 375a. u. et 
562. 5a. u.), although a small discontinuity is present at 
t = 375a. u. for the run B. 
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FIG. 5. Relative difference of dissociation probability for 
vs. time, between runs A (1 step), B (2 steps) and C (3 steps), 
| £-4 1 (blue) and |^| (red). 



One example of a relative difference, that for the tran- 
sition probability to state v = 13, is presented in FigJS] It 
varies inversely with the probability itself and it is always 
lower than 0.2% (except for negligible values of probabil- 
ity), despite the very small value of the final probability 
(4.47 x 1(T 6 ). 
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FIG. 7. Comparison of the relative difference on final dis- 



sociation probability (right bottom) \P 



B,C 



P£J/P£ SS , for 



t = 750a.u. and relative difference on transition probabili- 
ties to the 16 first bound states (other points) \P B ' C (\v)) — 
P A (\ V ))\/P A (\ V )), between runs A (1 step), B (2 steps) and 
C (3 steps), | B ~ A A | (green triangles) and | Z A \ (red points). 



Finally FigJT] shows a comparison of the final rela- 
tive differences (t = 750a. u.) of dissociation and transi- 
tion probabilities obtained with runs A, B and C. The 
precision is stable in relation to scale changes and it 
stays lower than 0, 2% once the probability in question 
is greater than 10~ 7 . Relative differences at very small 
probabilities are almost of the same magnitude as those 
at large probabilities. 



V. CONCLUSION 

We find that it is possible to drive the CATM with 
multi-step propagation, which allows us to decrease the 
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dimension of the Fourier basis used to describe the time 
variations. The multi-step propagation with the CATM 
is made possible by the development of a general form 
for the time-dependent absorbing potential, which can 
now constrain the Floquet state connected to any general 
initial superposition of states. The operating condition 
remains the same as for the simple case of an initial eigen- 
state of the free system, i.e. the time integral of the imag- 
inary part of the absorbing potential must be negative 
and large enough. The multi-step CATM preserves the 
two principal features of the one-step version: contrary to 
the standard wave packet propagations, it calculates the 
probabilities with an almost constant relative accuracy, 
whether large or small and it produces a dilatation of the 
i?F-spectrum in the complex plane and thus facilitates 
recursive treatment, by the RDWA of wave operator the- 



tential presented here is efficient enough to constrain any 
statevector to be transformed progressively into a given 
different one, as driven by the TDSE. This was the nec- 
essary condition to be able to propagate such an initial 
state with the CATM. The technical advantage of this is 
that it transforms a dynamical integration problem into 
an eigenvalue problem. 
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